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ABSTRACT 

We present results from two ~30 ks Chandra observations of the hot atmospheres of the merging 
galaxy groups centered around NGC 7618 and UGC 12491. Our images show the presence of arc-like 
sloshing cold fronts wrapped around each group center and ~ 100 kpc long spiral tails in both groups. 
Most interestingly, the cold fronts are highly distorted in both groups, exhibiting 'wings' along the 
fronts. These features resemble the structures predicted from non-viscous hydrodynamic simulations 
of gas sloshing, where Kelvin-Helmholtz instabilities (KHIs) distort the cold fronts. This is in contrast 
to the structure seen in many other sloshing and merger cold fronts, which are smooth and featureless 
at the current observational resolution. Both magnetic fields and viscosity have been invoked to 
explain the absence of KHIs in these smooth cold fronts, but the NGC 7618/UGC 12491 pair are two 
in a growing number of both sloshing and merger cold fronts that appear distorted. Magnetic fields 
and/or viscosity may be able to suppress the growth of KHIs at the cold fronts in some clusters and 
groups, but clearly not in all. We propose that the presence or absence of KHI-distortions in cold 
fronts can be used as a measure of the effective viscosity and/or magnetic field strengths in the ICM. 



1. INTRODUCTION 

Minor mergers are a common phenomenon in the 
growth of galaxy groups and clusters and leave charac- 
teristic observable features in the intra-cluster medium 
(ICM). The passage of a less massive subcluster through 
the main cluster offsets the central gas peak from the 
central potential well. Subsequently, the offset ICM falls 
back towards the potential minimum and thus oscillates 
- or sloshes - inside the cluster cent er ( Markevitch et al. 
2001 Ascasibar fc Markevitch|2006 ). The transfer ot an 



gular momentum in oil-center mergers gives the slosh- 
ing gas a spiral-like appearance. The sloshing forms 
cold fronts (CFs) in the gas, i.e. arc-shaped disconti- 
nuities in density, temperature, metallicity and thus X- 
ray brightness wrapped around cluster cores. Similar 
features can occur in equal-mass mergers if the merger 
partners pass each other at a sufficiently large distance 
such that their cores are only disturbed but not de- 
stroyed. Sloshing CFs are commonly seen i n clusters 



( Markevitch et al. 2003 Ghizzardi et al.|[2010[ ) and have 



been observed in detail in several clusters and groups (see 
also r eview by Markevitch & Vikhlinin 2 007|), e.g. Abel l 
1795 (IMarkevitch et al.|200l[|Bourdm fc Mazzotta |2008[), 
Abell 2142 jMa rkcvitch jF1i.||2000| |Owers_etjljl200^ 
Abell 496 ( |Dupke eUal.|2007|), N G(j j098 jp^all et al.| 
20091, NGC 5044 JUastamio et al.| | 2009l, NGCSSSff 
(TCTa chacek eTal] |2010[ ), NGC 5846 lM^h^ek_e£aL 
20 11 C j the P erseus cluster ( Churazov et al. 2003 ; S anders 
et al.||2005|), and the Virgo c l uster (Kimionescu et al. 



2010| ). ^Roediger et al.| ( |2011| |2012[ ) performed hydro 



dynamic merger simulations specifically tailored to the 
Virgo cluster and Abell 496. By varying the merger ge- 
ometry and mass ratio, the observed positions of the CFs 
and contrasts across them can be reproduced quantita- 
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tively and thus the recent merger history be constrained. 

Sloshing CFs have generally been considered as smooth 
arcs (e.g. in Virgo, Abell 2142). This is in contrast to pre- 
dictions from high-resolution purely hydrodynamic simu- 
lations, where the CFs are distorted by Kelvin-Helmholtz 
instabilities (KHIs) due to shear flows along th e front 



( ZuHone et al. 2010 Roediger et al. 2011 2012 1. Both 



viscosity and magnetic fields have been i nvoked to ex- 



plain the absence of KHIs in observed CFs ( |Keshet et al.| 
20101 jZuHone et al.|[20lT[) However, some CFs do have 
significant substructure. For example, the multiple slosh- 



mple, the multipL 
ing CFs in A496 clearly display a remarkable boxy mor- 
phology which had been unexplained unt il recently. Our 
high- resolution non- viscous simulations (Roediger et al. 
|2012[ ) reproduce these boxy features and demonstrate 
that they arise due to KHIs at the CFs plus their in- 
terplay wMi_j}rojej3t^^ CF in RX 



J1720 ( |Markevitch fc Vikhlinin|20 07) is somewhat boxy- 
shaped, and the CFs in NGC 7618 and UGC 12491 pre- 
sented in this paper are clearly distorted. 

NGC 7618 and UGC 12491 are the central galax- 
ies of two nearby, roughly equal mass galaxy groups 
(z=0.017309, d L =74.1 Mpc, 1' = 21 kpc). An 



ASCA/GIS image of the pair ( |Kraft et aL||2006[ ) shows 
extended X-ray emission centered on each of the domi- 
nant early- t ype galaxies. A s hort Chandra observation of 
NGC 7618 ( jKraft et al.|2006] found a CF - 20 kpc north 
of the nucleus and a spiral-like tail. We obtained new 30 
ks Chandra/ ACIS-S images of the central regions of both 
groups that confirm the CF plus spiral tail structure in 
NGC 7618 and find a similar morphology in UGC 12491. 
These CFs, the clear spiral-shaped gas asymmetries or 
tails in bot h groups and their ne arly identical recessional 
velocities ( Huchra et al.|[l999[ ) suggest a group/group 
merger taking place almost entirely in the plane of the 
sky. This favorable viewing geometry and the proximity 



2 



Roediger et al. 2012 




TABLE 1 

Number of counts in elliptical regions shown in Fig. [4] 



Fig. 1. — Two co-added, background-subtracted and exposure 
corrected 30 ks Chandra/ACIS-S images of the NGC 7618 and 
UGC 12491 groups in the 0.5-2.0 keV band, smoothed with a 4 
arcsec Gaussian kernel. 



of the groups make this an ideal target to study the gas 
dynamics of the merger process. 

In this letter, we discuss the distorted morphologies 
of the CFs present in both galaxies and the KHI as the 
likely origin of the distortions. A detailed analysis of the 
dynamic and thermodynamic state of the gas in this pair 
will be presented in a separate publication (M. Machacek 
et al., in preparation). This paper is organized as follows. 
Section 2 contains a short description of our data. In 
section 3, we discuss the merger history of this pair of 
groups. Finally, in section 4, we present the distorted 
morphology of the CFs and discuss the implications of 
our results. 

2. DATA 

NGC 7618 and UGC 12491 were observed for 34.5 and 
33.2 ks with Chandra/ACIS-S (OBSIDs 7895 and 7896), 
respectively, as part of the HRC Guaranteed Time Ob- 
servation program (PI: S. S. Murray). The data were fil- 
tered for periods of high background. Intervals where the 
rate in the 10-12 keV band was more than 3<r above the 
mean rate were removed. After filtering, 32.15 and 31.1 
ks of good time remained for NGC 7618 and UGC 12491, 
respectively. Bad pixels, hot columns, columns along 
node boundaries and point sources were also removed. 
All data were processed with CIAO version 4.3. We used 
blank sky background sets appropriate for the observa- 
tion date and instrument configuration renormalized to 
agree with the observations in the 10 — 12keV energy 
band, where particle background is expected to domi- 
nate. For more detail on the reduction and analysis pro- 
cedures, see Machacek et al., in preparation. 

3. THE MERGER BETWEEN THE TWO GROUPS 

A mosaic of the NGC 7618 and UGC 12491 field is 
shown in Figure [l] Bo th groups show char acteristic 
gas sloshing features (see Roediger et al. 2012 for a de- 
tailed list), i.e. the arc-shaped CFs around the cores in 
a spiral-like fashion, and the spiral-shaped brightness 
excesses (or spiral tails) extending to several CF radii. 
NGC 7618 shows a clear arc-shaped surface brightness 
discontinuity 20 kpc to the north of the nucleus, and a 
spiral tail that curves from the south through the east 
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and north. UGC 12491 exhibits a very similar structure 
but at a rotated orientation - a surface brightness dis- 
continuity to the south-west of the nucleus and a spiral 
tail curving from east over north then west. Tempera- 
ture measurements using the mean energy in the Fe L 
peak (0.7 S, E < 1.2 keV) as a proxy of temperature 
( David et al.| [2009) confirm the brighter side of the sur- 
face brightness discontinuities to be the cooler one and 
thus the discontinuities to be CFs. The morphological 
comparison of the orientation of the CFs and the spi 



ral ta ils to sloshing simulations (Roediger et al. 2011 
2012 1 suggests that both galaxy groups passed each other 



such that UGC 12491 came from the east-south-east and 
passed south-west of NGC 7618. 

We estimate the time since the closest passage between 
the groups in two ways. First, the almost identical red- 
shifts of both group centers indicate a merger in the 
plane of the sky. The average group gas temperatures are 
~ 1.2 keV, corresponding to a sound speed of 560 km s -1 . 
Around pericenter passage in cluster mergers, the rela- 
tive velocity between both clusters varies between Mach 
1 and 2 (e.g. |ZuHone et al.||2010l |Roediger et al.||2011[) , 
merger partners on orbits with close encounters reach 
up to Mach 3 at pericenter passage. In the case of the 
NGC 7618/UGC 12491 merger, both groups are of sim- 
ilar mass, but their cores have not been destroyed but 
only disturbed, which argues against a close encounter. 
Hence, assuming an average relative velocity between 
Mach 1.5 (840 kms" 1 ) and Mach 2 (1120 kms^ 1 ) since 
their closest approach, the projected distance of 300 kpc 
translates into an age of around 0.3 Gyr. Second, as 
the outwards motion of sloshing CFs is mainly governed 
by the host potential, the distance of the CFs from the 
center depends mainly on the time since the closest en- 
counter. Although the velocity of the CFs can in prin- 



ciple differ from cluster to cluster, Roediger et al.| ( 2011 
2012 1 found similar CF velocities tor Virgo and Abel 



Only in more mass ive clusters, like A2029, do the 



3m: 

CFs progress out faster (Roediger & ZuHone 2012). In 
all simulated systems, the CFs move with approximately 
constant velocity. Virgo is the least massive of these 
three systems. Assuming that the CFs in NGC 7618 
and UGC 12491 move at about t he same velocity as in 
Virgo (55 kpc/ Gyr, see Fig. 7 in Roediger et al. 2011), 
the group-centric distance of 20 kpc of the CFs in both 
groups translates into an age of about 0.36 Gyr, consis- 
tent with the first estimate. 

4. SUBSTRUCTURE OF COLD FRONTS AND 
IMPLICATIONS FOR ICM PROPERTIES 

The most interesting features of NGC 7618 and 
UGC 12491 are the irregular shapes of their sloshing CFs. 
In Figuresj2]and[3]we show the central regions for both 
galaxies where the color scale has been chosen to high- 
light the substructure along their CFs. The CF north of 
the nucleus of NGC 7618 displays a "nose" and a wing as 
labelled in Fig. [2j There may be more distortions/ wings 
along the western edge of the tail, south of the structure 
labelled "wing". UGC 12491 has a smooth CF to the 
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Fig. 2.— Chandra/ACIS-S image of NGC 7618 in the 
0.5-2.0 keV band, background-subtracted, exposure corrected, 
Gaussian-smoothed to 6 arcsec. The logarithmic color scale is cho- 
sen to highlight the substructure of the cold front. Prominent 
features are labelled. The dashed arc marks the cold front. 
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Fig. 3. — Same as Fig. [2] but for UGC 12491, Gaussian-smoothed 
to 4 arcsec. 

south-west which, however, terminates in two wings, one 
to the south-south-east, one to the north- north- west. Its 
tail appears to be split at about 30 kpc north-east from 
the nucleus and possibly again at 50 kpc to the north. 
Alternatively, these splits in the tail could be regarded 
as wings or distortions along the outside edge of the tail. 
In both groups, these CF substructures have linear scales 
of about 15 kpc. 

To demonstrate significance of the nose feature at 
NGC 7618, we compared the brightness of this feature 
with the adjacent background level. We performed this 
analysis on the raw image in the 0.7 to 1.4 keV band, 
where the ratio of source to background counts has been 
optimized only by the choice of energy band. We placed 
the elliptical region 1 shown in Fig. |4]over the nose. El- 
lipses 2 and 3 are at the same distance from the galaxy 
center, and ellipses 4 to 7 are placed north of the CF. All 
ellipses have the same size. We ensured that all ellipses 
are not contaminated by the brighter emission inside the 
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Fig. 4.— Raw Chandra/ACIS-S image of NGC 7618 in the 0.7- 
1.4 keV band, Gaussian-smoothed to 6 arcsec. The ellipse 1 covers 
the "nose" feature, ellipses 2 and 3 are at the same distance from 
the galaxy center, and ellipses 4 to 7 are placed such that they 
cover patches of high background. All ellipses have the same size, 
the number of counts in them is listed in Table ^ 

CF. The number of counts in each ellipse is listed in Ta- 
ble [T] Ellipse 2 is the brightest of the background ellipses 
ancTcontains 42 counts. Assuming Poisson noise implies 
a standard deviation of a = 6.5. The 60 counts in ellipse 
1 on the nose feature are 2.7 a above the background 
level as defined by ellipse 2. The average background 
level of regions 2 to 7 is 35 counts, corresponding to 
a = 5.9. Therefore the counts in the nose regions are 
4 a above this level. The length scale of the confidently 
detected structures is about 15 kpc. With the current 
data, smaller structures cannot be detected at compara- 
ble significance, because, assuming a comparable surface 
brightness also in smaller structures, the significance, 
i.e. the excess of source counts above the background 
divided by the standard deviation of the background, is 
proportional to the considered length scale. To better 
define the CF shape, deeper observations are required. 

We suggest that the distortions in these CFs are the 
result of KHIs, which are expected to arise due to shear 
flows along the CFs and are routinely seen in n on-viscous 
hydrodynamical simulations (e.g. Fig. Al in |Roediger 
et al.l pOlT] also |Roediger etal] [2TJ12{ |ZuHone et "aT 
2010p . ^soutlined in Sect. [I] ab ove, both distorted CFs 



and smooth, arc-like ones have been observed, and the 
presence or absence of KHI-like distortions constrains 
the effective viscosity of the ICM and tangential mag- 
netic fields along the fronts, which can both suppress the 
growth of the KHI. After a short remark on the effect of 
gravity we discuss these two IC M properties below. 

Gravity suppresses the KHI ( Chandrasekhar 1961) at 
wavelengths above length scales of 
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where we inserted typical values for the density contrast 
(pi and pi are the densities of the gas interior and exte- 
rior to the CF, respectively) across the CF, D, and the 
shear velocity U, along the front. We derived the gravita- 
tional acceleration g from the azimuthally averaged ICM 
density and temperatur e profiles in UGC 12491 given in 



the ACCEPT sample (Cavagnolo et al. 

This iengt. 



2009} , assum- 
scale A max is 



ing hydrostatic equilibrium, 
comparable to the cluster-centric distance of the CFs, 
i.e. the size of the low entropy core, which already limits 
the maximum perturbation length scale. Consequently, 
gravity does not play an important role in suppressing 
KHIs because the size of the cores is of the order of or 
smaller than A max . 

The KHI can be suppressed by sufficiently strong mag- 
netic fields aligned with the interface. The presence of 
KHI puts an upper limit on the magnetic field strength 
in these groups, because, if the combined magnetic pres- 
sures of the tangential magnetic fields at the hot and 
cold side of the CF, (B^ A +B% „,.)/(8ir) exceeds the limit 



(Vikhlinin & Markevitch 2002) 
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the growth of the KHI should be suppressed (picm is the 
ICM pressure, M is the Mach number of the shear flow, 
and Thot and T co \& are the temperatures at the warmer 
and colder side of the discontinuity, respectively). For 
an ICM pressure of picm = 10 -3 keVcm -3 at the CF 
this corresponds to a total magnetic field of 2 fiG at the 
CF. The presence of fields just weaker than this limit in- 
creases t he growth time of the instability, but only mod- 
erately. |ZuHone et al. ( 2011[ ) showed that gas sloshing 
typically amplifies the magnetic fields at the CFs by up 
to an order of magnitude, implying initial field strengths 
of around ~0.3 pG in these two groups, which is within 
other observational limits for groups (e.g. Guidetti et al. 
20l0|). 



Junk et al.| ( 2010 1 consider the KHI in the presence of 
viscosity and derive the dispersion relation for a shear 
flow of relative velocity, U, between two layers of incom- 
pressible fluids with constant densities, p\ and p^ with 
Pi > P2- They assume a constant kinematic viscosity, v, 
in both fluids. The corresponding growth timescale as a 
function of perturbation scale, A, can be expressed as 



1 



-(A) = -r [Vl + (A/L) 2 + l] with 



L-- 
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. TTV 
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The presence of viscosity separates two regimes of scale 
lengths above and below L. For A L the viscosity is 
irrelevant and the instability grows at the same rate as 
in the inviscid case. For A < L, the growth timescale is 
approximately constant at t w To, whereas in the inviscid 
case the growth timescale scales linearly with A. Thus, 
the viscosity slows down the growth rate compared to 
the inviscid case for A < L. Using the definition of L in 



Eqn. [5] this inequality can be written as 

XU 
v 



< ttVA 



(6) 



where ^ is the Reynolds number associated with the 
perturbation length A. Naturally, the viscosity becomes 
relevant at small Reynolds numbers. In fact, the constant 
growth time for A < L derived above is true for the 
onset of the instability only. Numerical tests for the long- 
term evolution of the KHI in viscous gases show that 
the viscosity quickly erases the shear flow parallel to the 
interface and thus slows down and eventually prevents 
the growth of the instability at small Reynolds numbers 
(Roediger et al., in prep.). 

We calculate the numerical values for L and To for CFs 
in different clusters in Table [2] We derive the kinematic 
viscosity v = M from the Spitzer value for the viscos- 



ity, p ( Spitzer |1956 ), assuming a Coulomb logarithm of 
In A = 40, ICM densities and temperatures as measured 
at each CF, and shear velocities either derived from sim- 
ulations (for Virgo, A496) or assuming shear velocities 
of Mach 0.4, a typical value found in sloshing simula- 
tions. Additionally, we list for each CF its distance to 
the cluster center and its age as derived by simulations 
(for Virgo and A496) or from the cluster-centric distance 
and a constant CF velocity as described above. Abell 
2142 is so mewhat hotter than t he ma ssive cluster sim- 
ulated by Roediger & ZuHone (2012). Assuming that 
the CFs in Abell 2142 move outwards with the same ve- 
locity as in this simulated cluster, we estimate that the 
sloshing in Abell 2142 was triggered about 2 Gyr ago. 
Sloshing CFs generally span more than 90 degrees in po- 
sition angle and thus more than their tqf in azimuthal 
extent. 

The estimates of L and To (see Table [2]) imply that vis- 
cosity does not significantly suppress the growth of KHIs 
in groups and poor clusters, because it can affect only 
perturbations at spatial scales significantly smaller than 
the CF radius or azimuthal extent. Also the estimated 
growth times are much shorter than the age of the CFs. 
Consequently, we should observe KHIs in NGC 7618, 
UGC 12491 and A496, which all show distorted fronts. 
In NGC 7618 and UGC 12491, the KHIs may grow even 
more rapidly, because in this roughly equal mass merger 
the shear flows along the fronts may be stronger than 
the Mach ~0.4 assumed here, which are typical for mi- 
nor mergers. The northern front in Virgo should also 
exhibit KHIs, which may be hidden in the low spatial 
resol ution of the current shall ow XMM-Newton observa- 
tion ( |Simionescu et aL 2010). The estimated values of 
L and To tor the hot cluster Abell 2142 are largest both 
in absolute terms and in comparison to the CF radii, 
and the viscosity is expected to have the strongest effect 
here. Indeed, there is no obvious evidence for KHI in the 
current Chandra data. 

A similarly diverse picture arises regarding the struc- 
ture or smoothness of merger CFs, i.e. the interfaces be- 
tween the atmosph eres of a cluster and an infalling sub- 



cluster or galaxy. Vikhlinin fc Markevitch| ( 2002 1 find 



that a tangential magnetic field ot about 10 pG is re- 
quired to_2reyen^l^ffl_a^he^mooth merger CF in Abell 



3667. Dursi & Pfrommer (2008) have shown that such 



magnetic draping layers can form. However, |Churazov 
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TABLE 2 

Full Spitzer viscosity should increase the growth time t(A) of the KHI for perturbations with length scales A < L, to 
t = ro (Eqn. H. This table lists values for the limiting length scale L and the resulting growth time to at CFs in 

DIFFERENT CLUSTERS, ALONG WITH ASSUMED VALUES FOR TEMPERATURE T, DENSITY 71 AND SHEAR VELOCITY U AT EACH CF. 

Furthermore, we list for each CF the distance to the cluster/group center rcF and its age tcf- 
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& Inogamov (20041 argue that the absence of the KHI 
near the stagnation point can be due to a small, un- 
observable intrinsic width of the curved fron t and the 
spatial varia tion of the shear flow along it. |Mazzotta| 
et al. ( 2002 ) suggest that the wings of the merger CF in 
A3bb7 indeed break up in KHIs. Recent deep Chandra 
and XMM observations of the merger fronts at the ellipti- 
cal ga laxies NGC 4552 jMachacek et al.|2006[ ) and NGC 
4472 ( Kraft et al.||2011[ ) tailing into the Virgo cluster re- 
veal significant substructure, i.e. horns, kinks and noses 
for both. In contrast, the upstream edge of NGC 1404 
( Machacek et al. |2005 ) falling into the Fornax cluster is 
not distorted. Like the sloshing CFs, some merger CFs 
are also distorted while others are not. The reason for 
this difference is yet unclear, it may be of dynamical na- 
ture, e.g. the amplitude of shear flow velocity or density 
contrast, or depend on the magnetic fields or transport 
coefficients of the ICM. 

Our simple estimate regarding the growth timescales 
of KHIs in the linear regime neglects several important 
issues that are relevant for real CFs: (i) The interface is 
not plane-parallel but curved, and the shear flow along it 
may vary with position and time. This is especially true 
for the merger CFs. In sloshing CFs, the interface moves 
away from the cluster center, stretching growing pertur- 
bations. This increases the growth time because the per- 
turbation length increases and because the density at the 
interface decreases, (ii) The growth time scales refer to 
the initial linear growth phase only. The evolution of 
CFs takes place on much longer timescales, where the 
KHI enters the non-linear regime, (iii) For a Spitzer-like 
viscosity, the kinematic viscosity is not constant across 
the interface, but the ratio of Reynolds numbers, between 
the hot and cold layer, scales with temperature contrast 
to the power of 3.5. (iv) The effect of tangential mag- 
netic fields and viscosity may combine in suppressing the 
KHI. 

The recent and upcoming deep and high resolution ob- 



servations of merging and sloshing cold fronts provide us 
with a wealth of information regarding the presence or 
absence of KHIs at both sloshing and merger CFs and 
hence the amplitude of magnetic fields and viscosity in 
the ICM. Understanding the data requires a systematic 
study of group and cluster CFs to determine which con- 
ditions lead to the growth or suppression of the KHI. 
The list of complications given above makes clear that 
simple 1 or 2D analogues will not be sufficient to pre- 
cisely interpret the observations. Additionally, both the 
required magnetic field strengths and viscosity for the 
suppression of the KHI depend also on the dynamical 
conditions like local density and shear velocity. Simu- 
lations tailored specifically to the object in question are 
required to correctly interpret the data. The direct com- 
parison of the overall properties of the CFs, i.e. positions 
and contrasts across them, allows the merger history and 
thus the current and recent dynamical state of the ICM 
to be constrained, e.g. densities, temperatures and shear 
flow histories at the CFs. Consequently, these dynami- 
cal properties can be disentangled from the ICM viscosity 
and magnetic field strengths, and these can be measured. 
We will follow this strategy in a series of forthcoming pa- 
pers. 
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